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Abstract 

Zitterbewegung, as it was originally described by Schrodinger, is an nnphysical, non-observable 
effect. We verify whether the effect can be observed in non-inertial reference frames/curved space- 
times, where the ambiguity in defining particle states results in a mixing of positive and negative 
freqnency modes. We explicitly demonstrate that snch a mixing is in fact necessary to obtain the 
correct classical valne for a particle’s velocity in a uniformly accelerated reference frame, whereas 
in cosmological spacetime a particle does indeed exhibit Zitterbewegnng. 


1 Introduction 


Zitterbewegung (“trembling motion”) of a free relativistic electron described by the Dirac 
equation was theorised by Schrodinger in 1930 [1]. It has been subsequently understood that 
Zitterbewegung is actually an unphysical, non-observable effect. In one-particle (quantum 
mechanical) interpretation one can demonstrate this by performing specihc unitary trans¬ 
formations suggested by Foldy and Wouthuysen [2], The effect disappears in the Foldy- 
Wouthuysen representation, and this suggests that a naive superposition of negative and 
positive energy solutions in the original Dirac representation is erroneous. A full resolution 
comes from the consistent second-quantized treatment according to which any real physi¬ 
cal particle (antiparticle) carries only positive energy and thus no superposition of negative 
and positive energy states is possible. The effect of virtual particles is taken into account 
to describe the physical electron (mass, charge and wave function renormalizations) and, 
obviously, does not affect its free motion^. 

In this paper we demonstrate that Zitterbewegung may show up as a physical effect in 
certain reference frames which cannot be dehned globally on an entire spacetime. In such 
spacetimes, the dehnition of particle states become observer-dependent and thus ambigu¬ 
ous. Well-known examples are reference frames related to stationary observers in black hole 
spacetimes and uniformly accelerating (Rindler) observers in flat Minkowski spacetime. Such 
observers perceive the Minkowski vacuum state as being a thermal radiation coming from 
the event horizon [5, 6], which is due to the non-trivial mixing of positive and negative fre¬ 
quency modes of quantized helds. Similar mixing and the related particle production occurs 
in time-dependent (cosmological) spacetimes [7]. We demonstrate that precisely this mix¬ 
ing of negative and positive frequency modes causes Zitterbewegung of a Dirac fermion in 
cosmological spacetime, whereas in Rindler spacetime it leads to a one-particle expectation 
value for velocity which exactly reproduces its classical counterpart. In what follows we 
restrict our discussion to (l-l-l)-dimensional toy models, since this case is technically simpler 
and captures the essential physics behind the phenomenon. We consider other physically 
interesting spaces elsewhere. 

The paper is organised as follows. In the next section we set up our notations and 
show that no Zitterbewegung appears in inertial reference frames. Section 3 is devoted to 
the calculation of one-particle velocity expectation value for uniformly accelerated observer. 
Finally, in Section 4 we demonstrate explicitly the physical Zitterbewegung within a toy 
cosmological model. The last section is devoted to conclusions and some useful formulae are 
collected in the Appendix. 

^Despite this clear-cut picture some researchers still advocate Zitterbewegung as a real effect and some 
even claim its experimental discovery [3]. In our opinion these claims are erroneous (see also the criticism 
in [4]). 
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2 No Zitterbewegung for inertial observers 

For the sake of completeness and to set up our notations we start our discussion by demon¬ 
strating the absence of Zitterbewegung in inertial reference frames, (l-l-l)-dimensional 
Minkowski spacetime is parameterized by Cartesian coordinates [t, x) and the interval reads: 

ds'^ = dt^ — dx'^ . ( 1 ) 

The Dirac spinor in (1+1) dimensions^ is a two-component complex held ■0, which satishes 
the Dirac equation 

- 'm)ip = 0 , ( 2 ) 

where the (1 + l)-dimensional y-matrices satisfy the anticommutation relation 

{7^7'^} = 2^7^^l2x2. 


For this calculation, we will use the following representation for the y-matrices: 


/ = 


0 1 
1 0 


= 


0 1 

-1 0 


The Dirac equation is obtained as Euler-Lagrange equation from the Lagrangian: 

C = i'07^5^'0 — rmljilj. 


( 3 ) 

( 4 ) 


The canonical formalism leads to the Hamiltonian density 

l-L = 'Kdo'ijj — C = '0(—iy^c?! + m)'il) = •^yoV’, (5) 

where 7r{x) = and in the last equality was used the equation of motion (2). The 

same Hamiltonian is of course obtained from Noether’s theorem, in view of the translational 
invariance of the Lagrangian (4), which leads to the conserved energy-momentum tensor 

( 6 ) 


whose component is the Hamiltonian density (5). 

By solving Dirac’s equation (2) (see, for example. Ref. [9]), one obtains the normalized 

positive/negative frequency = \/p^ + solutions, which are the familiar plane 
waves: 


y,(±) _ ^ f y/^P 4- P ^ ^^:iujpt+ipx 


( 7 ) 


^For a detailed presentation of Clifford algebras and spinors in arbitrary spacetime dimensions, see, for 
example. Ref. [8]. 
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Thus, a free quantum field as described by an inertial observer in Minkowski coordinates 
reads: 

= J , ( 8 ) 

where the annihilation (creation) operators a(a^), b{b^), respectively for particles and an¬ 
tiparticles, satisfy the standard anticommutation relations: 

{a{p),a\p')} = S{p - p') , {h{p), h\p')} = 5{p - p') , (9) 


other anticommutators being zero. Particle states are dehned as excitations over the vacuum 
state |0m): 

a(p)|0M) = Kp)|0m) = 0 , Vp . (10) 

In particular, particle and antiparticle with momentum p are described respectively by the 
states, 

\p) = aHp)|0M) , Ip) = &Hp)|0m) , (11) 

which are orthogonal to each other, (p|p) = 0, and both carry positive energy E = +\/p^ + 'm?. 

Next, we dehne a particle’s velocity as follows. The Dirac Lagrangian (2) is invariant 
under U{1) global transformations, which leads by Noether’s theorem to the conservation 
of the current Upon quantization, the current becomes an operator whose 

components are the current and charge density operators: 

^ ( 12 ) 

f = ; , (13) 


where : ... : denotes normal ordering of operators and p = {0,1} are the Lorentz vector 
indices in (l-l-l)-dimensional spacetime. Expressing the operators in terms of creation and 
annihilation operator, it is then easy to show that the velocity of a particle with momentum 
p and energy Up is: 


{p\j^\p) ^ 
{p\f\p) 


(14) 


This expression is just the classical velocity of a relativistic particle. Hence, no Zitterbewe- 
gung term appears within the consistent description of the Dirac particle in inertial reference 
frames. 


3 No Zitterbewegung in an uniformly accelerated ref¬ 
erence frame 

In this section we consider Dirac fermions in a uniformly accelerated reference frame. An 
observer moving with a constant positive acceleration, a, along the x-direction in (1-1-1)- 
dimensional flat Minkowski spacetime, accounts for an event horizon. Such a non-inertial 
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frame is described by the Rindler spacetime, with the following line element: 


ds‘^ = a^x^dr'^ — dy^ . 


( 15 ) 


For details about quantum held theory in Rindler space, see, for example. Ref. [10]. The 
Rindler coordinates (r, x) can be written in terms of Cartesian coordinates {t,x), which 
parametrize the hat Minkowski spacetime, as follows: 


T = — arctanh 
a 



, X = V. 


(16) 


It is obvious from the above equations that the Rindler coordinates cover only the |x| > t 
quadrant of the full Minkowski spacetime and x ^ [0,+cx3]. 

The Dirac equation in Rindler spacetime reads: 


i'y^dr = —iaxTldy^i’ ~ 


la 


7 ° t ‘ - axmi’ | iji 


( 17 ) 


The positive/negative frequency solutions to the above equation can be found in terms of 
modihed Bessel function of the second kind: 




mcosh(^f2) 

2 * 




iK^in-Umx) 


(18) 


where we have adopted the normalisation condition: ^ dx = 5 {fl' — D) and 

D > 0. The mathematical details of these calculations can be found in the Appendix. 

Given the mode solutions (18) we construct the Rindler quantum held as: 

/ OO ^ 

dn , (19) 

where A{B) and Al(i?I) are particle (antiparticle) creation and annihilation operators, re¬ 
spectively, acting on the Fock space as dehned by the Rindler observer, i.e., A^IOr) = 
-BoIOr) = 0, VD, |0r) being the Rindler vacuum state, etc. We then dehne the Rindler 
current and charge density operators as: 

3-1 = ax : ■■ , (20) 

30 = ; , ( 21 ) 


where the normal ordering is dehned with respect to the Rindler operators. Note also that 
the time evaluation is assumed with respect to the Rindler coordinate time r, which is related 
to the proper time Tp as axdr = dxp. 
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According to (16), the uniformly accelerated frame can be viewed as a Lorentz boosted 
inertial reference frame with an instantaneous velocity ar. Hence we can relate the boosted 
Minkowski modes to the Rindler modes via Bogoliubov transformations, i.e.: 



poo 

/ dn an^p + (3n,p , 

Jo 

(22) 

e-3T'V“^^G)(t(r,x),x(r,x)) = 

poo 

/ dn a*^^pZn^ + /d^^pZZ^ . 

Jo 

(23) 


Using these equations and the relations given in the Appendix we have computed the Bo- 
goliubov coefficients (see the details in the Appendix): 


1 — i 




iClOp 


2a7Tm cosh(0p) cosh(-7r) 




In (24) we have parameterized the Minkowski frequency and momentum as: 

ujp = mcosh(0p) , p = msinh(dp) . 


(24) 


(25) 


Given the above Bogoliubov coefficients we can map Rindler (A,R) and Minkowski (a, b) 
operators: 


= JZc dk [an,A + , 


(26) 




= IZo dk {f3n,kak + an^Zl) ■ 


(27) 


As described in the Appendix, we hnally obtain the expectation values for the Rindler current 
and charge density operators given generically by (A.3.1) (see also for the specihc case of 
the current density (A.3.2)), respectively, for a one-particle Minkowski state \p): 


{p\jk\p) = 






' —CO J —OO 


axsinh(6*p — ar) 
TT cosh(6*p) 


(28) 






<^h,p0^ft',P 


cosh (6'p — ar) 
7rcosh(6*p) 


(29) 


Importantly, the limit of integration in the above equations reflect the fact that both the 
positive and negative frequency solutions do contribute to the hnal result. Also we note that 
the computed expectation values are hnite due to the Rindler normal ordering adopted in 
(20) and (21). With these in hand, we compute the velocity of a particle in a uniformly 
accelerated reference frame. 


vr 


{p\Jr\p) 

{p\Jr\p) 


tanh(0p — ar) , 


( 30 ) 


which is precisely the classical velocity. Hence, we see that, similar to inertial frames, there 
is no Zitterbewegung effect in uniformly accelerated reference frames. 
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4 Cosmological Zitterbewegung 

Let us consider now a toy (l+l)-dimensional cosmological model (see. e.g., Ref. [10]) with 
the line element 

= a^{ri){drf — dx"^), (31) 

which describes a change of the scale factor from a constant, Uq — y, in the infinite past 
{rj —>■ —cxd) to another constant ao + x in the infinite future (p —>■ cxo): 

a{r]) = ao + xtanh(£r 7 ). (32) 


This model describes an asymptotically static spacetime which undergoes a period of smooth 
expansion. The Dirac equation in this spacetime now reads 


2 a{r]) 


maij]) 


V’ = 0 , 


(33) 


where the prime denotes the derivative with respect to r]. In the following calculations we 
shall use the notations: 

^injout ~ (®0 T X)) 


^ in lout = \jv^ ^ 


2^2 

in!out ’ 




(34) 


The solutions to this equation can be found using the Dirac representation for the y-matrices: 
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0 


0 -1 y ' V -1 0 / 


(35) 


The positive energy solutions which match the flat-space plane waves in the infinite past, 
denoted below as “in” states, are found to be (see Appendix for a brief discussion): 




sh(e? 7 )) ®(d) 


\/2 




Ui 

U2 


^ipx 


(36) 


with 


U\/^in -^1; 

U 2 = — {!jJin - main)~^ 

p 


X 


Fi 

F 2 


—mao + tanh(e77)(ze — mx) + - [{i<^in — e)(tanh(£? 7 ) — 1) — {iuJout + e)(tanh(£r 7 ) -|- 1)] 


lUJi. 


[4z [(a;_)^ — m^x^] + 4£(a;_ — mx)] (tanh(£77) -|- 1) (tanh(£r7) — 1) 


( 37 ) 
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Here we have used the following condensed notation: 


= 2^1 - (w- - mx ), 1 + -(ca_ + mx); 1 
e e 


^tanh(£77) + , 

" ' 1 1 

F 2 = 2^1 ( 2 + - {u- + mx ), 1 + -(a;_ - mx); 2 - - tanh(£r/) + - 

e e e 2 2 


( 38 ) 


with 2 Fi{a,b; c; z) being the ordinary hypergeometric function. The states which solve the 
Dirac equation and match the plane wave solutions in the inhnite future are denoted as 
“out” states and the positive energy solutions are found to be: 


jIj+ ^ — ("gE*? _|_ g-£f(^g-*‘-^+»(g-(l+^)ln(2cosh(e?j)) 


a{r]) 2 fUi 


yJ^UJout \^2 


^ipx 


(39) 


with 

U\ \/UJout T F\ 

U2 ix^out ^TT'^out) ^ 
P 


X 


Fi 

h 


—mao + tanh(£77)(f£ — mx) + - — e)(tanh(£77) — 1) — {iuJout + e){tanh.{er]) + 1)] 


[4i — m^x^ + 4£(a;_ — mx) {t&nh.{er]) + 1) (tanh(£? 7 ) — 1) 


Here we have introduced the following notation 
-^1 = 2F1 


i r N , * / X , 1 1 , / N 

-{uj_ - mx), 1 + -(w- + mx); 1 +-; - - - tanh(£r/) 

e e e 2 2 


F2 — 2F1 


2 + - (n;_ + mx ), 1 + -(w- - mx); 2 + ^ ^ tanh(£77) 

£ £ e 2 2 


(40) 


(41) 


The negative energy states can be found by using the charge conjugation operator as follows 

. (42) 

In our 2-D Dirac representation, the charge conjugation operator is 


C = 


0 1 
1 0 


The held operator can be expanded in terms of the above mode functions as follows: 

^in = I + b\-p)i^l~X) > for r; < ?7o , 


^out 




-) 

out 


, for ?7 > 770 • 


(43) 

(44) 

(45) 
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The two sets of creation and annihilation operators in the above equations are related to 
each other via the Bogoliubov transformations: 


A{p) = a{p)a{p)+l3*{-p)h\-p) , 


(46) 


B\p) = l3{-p)a{-p) + a*{p)b\p) , 

where 

^ ^/^t^/ooin + maj _ r(l - i^)r(-^) 

+ maout r(l - -^{u+ - mx)r(-^(a;+ + mx)) 

- bl + maj _ r(i - _ 

P y/^in\/^out 'ITT'^^out r(~(ca_ ?77.^)r(l T —{oj— T ^X)) 

Assume now that we prepare a one-particle (Heisenberg) state with momentum p at some 
earlier time p —)■ —oo: \p) = a'*'(p) |0M;m), and we wish to compute the velocity of this state 
at some later time p ^ oo, that is 

(pIjS,.iIp> ’ 


(47) 

(48) 

(49) 


where 

jLt =■ ^lut^^Y^out ■ 

The velocity can be calculated generically, and is found to be: 

I cosh(^mx) 


IpI 

y — _iiJ_ 

^out 


ta,nh{jU!out) tanh{ju}i„) sinh(^a;o„t) sinh(^aji„) 


+ ^ + e 


(50) 


(51) 


where we have introduced 
sinh(^!^) 


X sinh (w+ -l- my) j sinh (ta+ — my) j sinh (w- -|- my) j sinh (ta_ — my) j 


xT r ( --(a;+ - mx) ] T ( --(a;+ + mx) ] T ( --{u_ + mx) ] T ( --{u_ - mx) 


V ^ J 


(52) 

We can consider the case where the spacetime grows rapidly at one instant in time. This 
case can be evaluated by taking the limit as e —)■ cx), where the spacetime approaches a step 
function and the in-states are instantaneously converted to out-states at time p = 0. In this 
limit, the velocity becomes: 


V = 


\p\ (lin(^out 
^out ^in^out 


' P'^ , 2 

-h m 

1+ ( cos {2ujoutp) 

_Ctin^out 

\ ^in / 


(53) 



As envisaged earlier in [11], we explicitly enconnter in (53) the time-dependent oscillatory 
term, which is very similar to Schrodinger’s Zitterbewegung term. The effect vanishes for 
massless particles as in the standard case. As a consistency check, we note that for = aout, 
the Zitterbewegung term disappears from (53) and one recovers the usual velocity formula 
given in (14). 

We can further consider the case when the Compton wavelength of the particle is much 
larger than the width of the spacetime change i.e — <^ 1. In this limit, maintaining only 
linear terms in y, we find that the Bogoliubov coefficients become 


Oip 

Pp 


■•»((?)')■ 

---ifntofF)’) - 

£ p smh(7r|) \ \ e / ) 


( 54 ) 


and hence for the velocity expectation value (keeping quadratic ^ terms for clarity), we 
achieve 


V = 


bl 

^out 


1 + 


m\ 2 


TT 


e J 2 sinh(—) 


(^(Xcnit^in ^out ^in) , ^outi^^out ^in)^ \ 

H- cos{2uoutV> 


2 sinh(—) 


P 


( 65 ) 


In this small mass limit, it is clear that the classical velocity is achieved and only 
oscillatory terms of order (y)^ contribute. 


5 Conclusion and outlook 

In this paper we have discussed the Zitterbewegung phenomenon in non-inertial reference 
frames /curved spacetimes, where the positive and negative frequency modes of a quantum 
field get mixed. We have explicitly demonstrated that this mixing is essential to obtain the 
classical velocity in the Rindler spacetime, that is, similar to an inertial observer, Zitterbe¬ 
wegung is an unphysical and thus unobservable effect for an uniformly accelerated observer. 
Within a toy cosmological model, however, we have found that Zitterbewegung can indeed 
occur. 

It certainly is interesting to study other non-inertial reference frames/curved spacetimes 
in which Zitterbewegung manifests itself as an observable phenomenon. Besides potential 
applications to cosmology, this study may offer an intriguing possibility to verify Zitterbewe¬ 
gung in laboratory experiments that explore analogue curved spacetimes (see, e.g., [12] for 
a recent review). The obvious challenge for such experiments is to maintain a low enough 
temperature, so that thermal fluctuations do not overshadow the quantum Zitterbewegung 
effect. In this regard, trapped ions may constitute a promising experimental set-up for 
the observation of the cosmological Zitterbewegung effect, providing rapid and controlled 
expansion/contraction of the ion trap can be achieved [13]. 
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A Appendix 


In this Appendix we collect useful formulas and relations used for obtaining the results 
discussed in the main text. 


A.l Normalisation of Rindler modes 

The solutions to the Dirac equation in Rindler spacetime (17) are: 


^(±) = I _ 


(A.1.1) 


We calculate the normalisation factor, Nq, by imposing the normalisation condition: 

/•rv~i 


Jo 

To evaluate the normalisation factor, we solve the following integral: 


dz K i{z)K-n_i{z) = 


in 


jQ a s as 4sinh(^(D'— D)) cosh(^(D'-|-D)) 

using the integral relation 6.576(4) for modified Bessel functions from [14]: 


dzzPKy{z)K^>{z) = 22-PTii+p) ^ 

l+P+('»+'ij') l+p-('iJ+'»0 1+P+('»-'»') l+p-{v-v ') ' 


along with the following Gamma function relations [14]: 

F(z*) = F*(^) , 

zF(z) = F(z 1) , 


|F(z.)p = 


n 


2;sinh(7r2;) ’ 


1/1 \ IO ^ 

ru + iz t = —rrv’ 

2 cosh(7r2;) 


|r(i + z:.)P = - 


nz 


sinh(7r2;) 


(A.1.2) 


(A.1.3) 


(A.1.4) 


(A.1.5) 
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The integral (A. 1.3) is one of two terms that arise from the normalisation condition (A. 1.2). 

The second term is obtained by taking —)■ fl'. The total contribution is identically 0, 

except at the pole = fl'. We can evaluate the pole by appealing to the Sokhotski-Plemelj 

formula. In a distributional sense it can be shown that: 

1 , P.V. 

lim — ---= —i2a6{Q — fl) H- , , _ , (A.1.6) 

£->o smh(^(fl' — fl)) + smh(^(f2' — fl)) 

where P.V. stands for the Cauchy Principal Value. By combining (A.1.1) - (A.1.3) with 
(A. 1.6) we find the normalisation constant: 


No, = 


m cosh(^ 7 r) 
TT^a 


(A.1.7) 


which gives us the mode solutions (18). 

A.2 Bogoliubov coefficients 

Here we present the explicit calculations only for the ao. coefficient, as the (io coefficient is 
obtained trivially by changing the sign of as in (24). Utilizing the normalisation condition 
on the mode functions we can rewrite (22) as 




d-x 


(A.2.1) 


Written out explicitly, this is: 


mcosh(—)\ / K_-n,i 

dx I -I I ^ 1 


7ra 


a 2 , 


(27r)2 y^2^ 




h{ar)-\-ipx cosh(ar) / x/Up-p 


We then use the change of variables: 

u = mcosh(0p) , k = msinh(0p) , 

to obtain: 


(A.2.2) 
(A.2.3) 


an,e,= I ^ a ) ^inr^-imxsinHar-e,) 


2n 


naup 


This can be solved using the integral 6.611(3) [14]: 

a + xj 


Ar_in + i(mx)-fA:_^n + i(mx)e 2^“^ 

a 2 a ' 2 

(A.2.4) 


After some algebra, we obtain: 




a + xj — (5“^ 


1 — i 


ttO 

e a e 2 a 


oinfip — 


2a'Km cosh(0p) cosh(^7r) 


(A.2.5) 


(A.2.6) 


11 
























A.3 Rindler current and charge density calculation 

We are interested in calcnlating the expectation values (lp|j^|lp), where is given in (20,21). 
Explicitly, we have: 


(lp\ : : |lp^ = 


^OO POO . . + 


dQ dfl' 


— OO J —OO 

POO POO 




(+) 

p' 






d£L dVt' 


/o Jo 


dk Un'’] ' a-n,e^a_n\9A0) > 




(A.3.1) 

where O is 1 for the charge density and 'j^'y^ax for the current density. The third term arises 
from the normal ordering of the Rindler operators and exactly cancels with the divergent 
second term (this can be shown by using the relationship between and /?n, along with the 
standard normalisation of Bogoliubov coefficients). The hrst term is the only contribution 
to this expectation value. For the current density, this integral is 


bOO POO 


' —OO ^ —OO 


dVtdVt' 

) 

Writing this out explicitly we have: 


(A.3.2) 


^OO POO 


dVt dfl 


' —OO J —OO 


47r3acosh(0p) 




(A.3.3) 


where we have used the shorthand notation i (mx), K±± = (my) and 

a' 2 a 2 

primes indicate that the order of the Bessel function is a function of hi'. We now use the 
relation 6.664(6) from [14]: 


Ka(x)K/^(x) = 2 / dycosh(y(Q; +/5 ))A’q,_/3(2x cosh(j/)), 

Jo 

as well as the Bessel recurrence relation 8.486(10) [14]: 

2a 


(A.3.4) 


Ka+l{x) - Ka-l{x) = - Kn{x), 

X 


(A.3.5) 


to show that: 


K++K'+_ - R+-<+ = 2^ 


/g amx cosh(j/) 

After applying the change of coordinates 


(A.3.6) 


n — - (11+ + r2_), 
fi' = i (n+ - n_), 


(A.3.7) 
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the total integral becomes: 


foo /‘oo poo 


’ —OO J —oo J 0 




The f2+ integral gives two delta functions, specihcally: 


r I'y , , e f-y , 


(A.3.8) 

(A.3.9) 


The integral is of the form: 

i fil- 

/ -fl_K±a (2mxcosh(y))zsin ( — (Op + ar 

J-oo a “ \ a ' 

which can be solved using 6.795(3) [14]: 

/ OO 

X sin (ax) 7742,(6) do: = 7r6sinh(a)e“^“®*^*'“^ 

-OO 

After solving both the r2_ and the 72+ integrals we obtain: 
a [6{y + i^) + S{y - i^)] 


(A.3.10) 


(A.3.11) 


dy- 


27r2 cosh (dp) 


7rxsinh(dp - a^)e-2"^xcosh(y)cosh(0,4+ar)_ (A.3.12) 


In order to take the integral over y, we deform the integral to the line z = ±ia+y to integrate 
the complex arguments of the delta function. The hnal result is 


axsinh(dp — ar) 

TT cosh (dp) 

Following an almost identical procedure, the charge density (29) can be found. 


(A.3.13) 


A.4 Mode Solutions for the Cosmological Spacetime 

The general solution to the equation of motion, given in (33) can be written in the form 


_i faUi{r]) + 13 Ui{ri)[uin -)■ -Uin] 


' (aU2iv)+(3U,{y)[ul 




^ipx 


with 


Ul = 2^1 


2 % 2 1 1 

-{u- - my), 1 + -{u- + my); (1 - -Uin)] (- tanh(e77) + -) 

e e e 2 2 




(A.4.1) 
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and 


P 


X 


Fi 

F 2 


lUJr 


—mao + tanh(e77)(ze — mx) + - [{ioJin — £)(tanh(e?7) — 1) — (iuJout + £)(tanh(£?7) + 1)] 
— [4z [oj'i — m^x^] + 4£(a;_ — mx)] (tanh(£77) + 1) {tanh.{eri) — 1) 


(A.4.2) 

In order to find the in/ont states, we match this general solntion to the plane wave solntions 
in the limits t] ^ —00 and r] ^ 00 respectively. The plane wave mode solntions to for a 
spinor field in the 2D Dirac representation are: 


t 


0 ( 77 )- 


IpI •-I ^ 


'-^1°^^ ^ dWjn/outVr p' ^J^injout ^O,in/out 

^ ( ■y yj^inlout — fnain/out\ iuJi^/autV- 




injout 


\f^ in/out'^ V '\J^injout T '^^in/out 


(A.4.3) 


ipx 


The in and ont states given in (36) and (39) approach these solntions in their respective 

limits. 


A.5 Cosmological Bogoliubov Coefficients 

The Bogolinbov coefficients can be obtained by determining how the in state evolves as 
7] —)■ cxo. To do this we nse the following hypergeometric fnnction relation[14]: 

F(a, /3; 7 ; z) ^ a + /? - 7 + 1; 1 - 2 ^) 

r(7-a)r(7-/3) 

+ (1 - + ^ F{x - a, 7 - /?; 7 - a - /? + 1 ; 1 - z), 

r(a)r(^) 

(A.5.1) 

which, when applied to the in state hypergeometric fnnction (in the first component), gives: 


2 F 1 


^(cn. - mx), 1 + ^(n;_ + my); fl - j ; Qtanh(£r 7 ) + ^ 


^ p ^_ iuJout'j 


r(l - l{uj+ - mx)T{-Hu+ + mx)) 


2 F 1 


C N N . 'i'^OUt /I 1 , / N 

-(a;_ - mx), 1 + -(w- + my); 1 +-; 7-7 tanh(£r 7 ) 

e e e \ 2 2 


X 2 F 1 


• In 2 cosh £77 


r(i- ^ ^ i^out ^ 


r(^(n;_ - my)r(l + ^{u- + my)) 


1 - -{uj+ - mx),--{uj+ + mx); 1 - ^ ^ tanh(e 7 ) 

s s s \ 2i 1 


(A.5.2) 
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which allows us to write 


(A.5.3) 


i^tn = (^pi^p,out + M-p,out 
with Up and (5p defined as in (48) and (49). 
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